Confined Quantum Time of Arrival for Vanishing Potential 
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We give full account of our recent report in [E.A. Galapon, R. Caballar, R. Bahague Phys. Rev. 
Let. 93 180406 (2004)] where it is shown that formulating the free quantum time of arrival problem 
in a segment of the real line suggests rephrasing the quantum time of arrival problem to finding 
a complete set of states that evolve to unitarily arrive at a given point at a definite time. For a 
spatially confined particle, here it is shown explicitly that the problem admits a solution in the form 
of an eigenvalue problem of a class of compact and self-adjoint time of arrival operators derived by 
a quantization of the classical time of arrival. The eigenfunctions of these operators are numerically 
demonstrated to unitarilly arrive at the origin at their respective eigenvalues. 
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I. INTRODUCTION 



The quantum free time of arrival problem has been the 
subject of numerous investigations in the past and in cur- 
renttimes JÍffippiÍàalÍÉIÏÍIllllÍlïllïlE 

El El H3, WÍ HI ÜÜ] The past has been mostly marked 
by stark pessimism fueled by Pauli's well-known theorem 
on the non-existence of self-adjoint time operators [l|, 
and Allcock's resigned conclusion that no ideal time of 
arrival distribution exists within the Standard framework 
of quantum mechanics Q. The current times, in contrast, 
is marked by optimism fueled by the generalization of 
quantum observables to include positive operator valued 
measures (POVM) [IlElElIïlElEllISlIlIllIll, 
and the advent of Kijowski's ideal quantum time of ar- 
rival distribution jfj. This optimism has been further 
strengthened by the realization that Kijowski's distribu- 
tion can be completely derived from a POVM arising 
from the quantization of the classical free time of arrival 
[Ï9j . And this optimism has been capped by the recog- 
nition that Kijowski's distribution has an opcrational 
meaning [2ÏJ |3^ US 03 EÜ- These recent significant 
developments have had led to the impression that the 
quantum free time of arrival problem has been solved, 
and only its proper generalization and experimental re- 
alization are left undone |3f3 |. 

While the above developments have unfolded, an unex- 
pected development has been in the silent offing. First, 
it was demonstrated by one of us that Pauli's theorem 
does not hold within the single Hilbert space formu- 
lation of quantum mechanics, and showed the consis- 
tency of self-adjoint and bounded time operators canoni- 
cally conjugate with a semibounded discrete Hamiltonian 
[371 Í39t |40| . This in turn has led to the realization that 
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the non-self-adjointness of the quantized free time of ar- 
rival (TOA) operator has nothing to do with the semi- 
boundedness of the Hamiltonian |37| . as Pauli's theorem 
would assert otherwise. This then led to the unexpected 
result that the non-self-adjointness of the TOA-operator 
can in fact be lifted by spatial confinement. Thus the 
concept of confined quantum time of arrival (CTOA) was 
introduced [38|. The CTOA-operators form a class of 
compact and self-adjoint operators canonically conjugate 
with their respective Hamiltonians in a non-dense sub- 
space of the system Hilbert space. Being compact, the 
CTOA-operators posses discrete spectrum, and a com- 
plete set of mutually orthogonal square integrable eigen- 
functions. However, the interpretation of the spectral 
properties of the CTOA-operators was not clear. Thus 
in |41| we addressed the issue of interpretation and it 
entailed rephrasing the quantum time of arrival problem 
in finding a complete set of states that unitarily arrive 
at some predetermined point. The eigenfunctions of the 
CTOA-operators are found to be states that evolve to 
unitarily arrive at the origin at their respective eigenval- 
ues. This result has repercussions beyond the quantum 
time of arrival problem which we will address elsewhere. 

In this paper we give full account of the confined quan- 
tum time of arrivals for vanishing potentials. In Section- 
[l·ljwe give a short review of the time of arrival operator 
in the real line. In Section- lïïïl we give the detailed for- 
mulation of the confined quantum time of arrival oper- 
ators. In Section- ITVl we investigate the structure of the 
conjugacy of the CTOA-operators with their Hamiltoni- 
ans. In Section-lvlwe study the symmetries of the CTOA 
operators and from them find the general properties of 
their eigenfunctions and eigenvalues. In Section- IVIl we 
solve analytically the CTOA-operator eigenvalue prob- 
lem. In Scction- IVIII we study numerically the dynami- 
cal properties of the eigenfunctions and provide unam- 
biguous interpretation for the spectral properties of the 
CTOA-operators. In Section- IYlIII we discuss the physical 
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repercussions of our current results. III. THE CONFINED TIME OF ARRIVAL 

OPERATORS 



II. TIME OF ARRIVAL OPERATOR IN THE 
REAL LINE 

The QTOA-problem is traditionally the problem of 
finding the time of arrival (TOA) distribution of a struc- 
tureless particle prepared in some initial state at a given 
point, say, at the origin. This operator is presumed to 
be the quantized classical-TOA in unbounded free space. 
That is, if a classical free particle, of mass y, in one dimen- 
sion at location q with momentum p, will arrive at the 
origin at the time T(q,p) = —yqp^ 1 , then the quantum 
TOA-distribution must be derivable from the quantiza- 
tion of T(q,p), from the operator 

T = -^(qp- 1 + p^q). (1) 

Formally the time of arrival operator T is canonically 
conjugate to the free Hamiltonian, H = (2fi)~ 1 p 2 , i.e. 
[H,T]=ih. 

Much of the study on equation-Q) has been done in 
momentum representation, in which it assumes the for- 
mal form 

_ ihfx í 1 2 d \ 
2 \p 2 pdp) 

As disccused in 0, T is a densely defined, unbounded 
operator in Tí p := L 2 (R, dp). And it has the following 
degenerate non-square integrable eigenfunctions 

where a = ±1. This set of eigenfunctions is complete, 

i-e,Ea J"-!L dí ^°f V)$? \p) = 5 (P~P')- However, they 
are nonorthogonal, i.e. 

j™Jp^?(jp)^(p) = \5 aa , ^(t-tO + ip^) . 

This nonorthogonality is a reflection of the non-self- 
adjointness of T in the real line. In fact T is maxi- 
mally symmetric so that it has no self-adjoint extension. 
The space-time properties of the TOA-operator eigen- 
functions are studied in |26[ . 

For a long time the non-self-adjointness of equation-QJ 
in the real line has been construed as a consequence of 
Pauli's theorem, in particular, from the semiboundedness 
of the Hamiltonian. However, as we will show in the 
following section, the non-self-adjointness has nothing to 
do with the semiboundedness of the Hamiltonian. 



A. The System 

Let the particle be confmed between two points with 
lcngth 21. If p 7^ and \q\ < l, the classical time of arrival 
at the origin (the hrst time of arrival, i.e. arrival without 
reflection from the boundaries) and the Hamiltonian are 
still given by T = —fiqp^ 1 and H = (2/i)~ 1 p 2 , respec- 
tively; moreover, T remains canonically conjugate with 
the Hamiltonian. Then equation-Q is still the totally 
symmetric quantized form of T even when the particle 
is confmed, and it likewise remains canonically conjugate 
with the Hamiltonian. 

To give meaning to T for the spatially confmed par- 
ticle, we attach the Hilbert space Tí = L 2 [— 1,1] to the 
system. The position operator is unique and is given by 
the bounded operator q, (qip)(q) = qip(q) for all <p(q) 
in TL. On the other hand, the momentum operator and 
the Hamiltonian are not unique, and have to be consid- 
ered carefully. Our choice is dictated by the assumption 
of closedness of the system and the requirement of con- 
sisteney with quantization: We assume the system to be 
conservative and we require that the evolution of the sys- 
tem be generated by a purely kinetic Hamiltonian. The 
former requires a self-adjoint Hamiltonian to ensure that 
time evolution is unitary. The later requires a self-adjoint 
momentum operator commuting with the Hamiltonian to 
ensure that the quantum Hamiltonian is the quantization 
of the purely kinetic Hamiltonian of the freely evolving 
classical particle between the boundaries. 

One of the possible Hamiltonians that can be assigned 
to the system is the textbook Hamiltonian where the do- 
main of the Hamiltonian operator is restricted to those 
vectors that vanish at the boundaries. But this Hamil- 
tonian, while self-adjoint, does not satisfy the second re- 
quirement: No self-adjoint momentum operator commut- 
ing with the Hamiltonian exists. The reason for this is 
that the eigenfunctions of self-adjoint momentum opera- 
tor in a bounded segment of the line are plane waves and 
none of these eigenfunctions vanishes at the boundaries 
(see below). This means that the momentum operator 
and the Hamiltonian do not have a common set of eigen- 
vectors. The Hamiltonian then can not be purely kinetic. 
For this reason, we abandon this Hamiltonian and con- 
sider another. 

Now for every 7 in the interval (—n/2, < n/2], there 
exists a self-adjoint momentum operator given by the op- 
erator p 7 = —ihd q whose domain consists of those vec- 
tors <j){g) in 7í with square integrable first derivatives, 
i.e. J \4>'(q)\ dq < 00, satisfying the boundary condi- 
tion 4>(—l) — e~ 2l7 0(7). With p 7 self-adjoint, the kinetic 
energy operator K 7 = t^P 7 is consequently self-adjoint. 
Thus the Hamiltonian is purely kinetic, 

H, = ^ (2) 
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the domain of which consists of all vectors 4>(q) in the do- 
main of the momentum operator p 7 such (p 7 0)((j) is still 
in the domain of p 7 . The momentum and the Hamilto- 
nian then commute and have the common set of plane 
wave eigenvectors 



1 



exp 



i (7 + kn) 



(3) 



where k = 0, ±1, ±2 • • • ; and their respective eigenvalues 
arep 7i fc = H(-f + kir)^ 1 and E ltk = h 2 + kir) 2 (2^il 2 )~ 1 . 

One may ask which of these infinitely many Hamilto- 
nians should we use in constructing the operators corre- 
sponding to the confined classical time of arrival. We will 
find below that we have to consider the whole window 
(— 7r,7r/2] in order to cover the entire symmetry of the 
classical time of arrival enumerated above in the quan- 
tum domain. Likcwise, we shall show that for every 7 
Equation-(m defines a self-adjoint operator, T 7 , satisfy- 
ing the canonical commutation relation with the Hamil- 
tonian H 7 in a closed subspace of the Hilbert space. And 
that the operators T 7 can be legitimately called as time 
of arrival operators; this follows from our study of the 
dynamics of the cigcnfunctions of the TOA-operators. 



B. Non-periodic Confined Time of Arrival 
Operators 

Now let us consider T for 7 7^ 0. Since q appears in 
first power in T, T is an operator if the inverse of p 7 ex- 
ists. Since zero is not an eigenvalue of p 7 , the inverse p" 1 
exists, and is in fact bounded and self-adjoint. Then it 
follows that, for every 7, T is a bounded, symmetric op- 
erator. Thus T is self-adjoint. For a given 7, we identify 
T with the operator 



T 7 = -^(Wj 2 



p 7 1 q), 



(4) 



derived from the formal operator T by replacing p with 
p 7 . We shall refer to T 7 as the non-periodic confined 
time of arrival (CTOA) operator for a given I7I < ir. 

In coordinate representation, T 7 assumes the form of 
a Fredholm integral operator 

i·l 



(2»( 9 ) = / {q\T^\q')y{q')dq'. 



(5) 



where the kernel is given by 



/ 1 -r 1 '\ (g + gO / '. 



m(q - q') + e- l m(q' - q)) , 

(6) 

in which IL(q — q') is the Heaviside function (see 
Appendix- A for the derivation of the kernel) . The kernel 
(q\ T 7 \q') is both symmetric and square integrable, i.e. 
(q\ T 7 \q') = (q'\ T 7 \q)* and J J \(q\ T 7 |ç')| 2 dqdq' < oc, 
respectively. This reaffirms the self-adjointness of T 7 . 
These two properties imply that T 7 is a compact, self- 
adjoint operator: It possesses a complete set of eigenfunc- 
tions with a corresponding discrete set of eigenvalues. 



C. Periòdic Confined Time of Arrival Operator 

For the periòdic case 7 = 0, p _1 is ill-dcfincd because 
p has no inverse, the zero being an eigenvalue of p. But 
this can be remedied. The pathology arises from the one- 
dimensional subspace spanned by the state of vanishing 
momentum, the null subspace A/*(p). But this subspace 
has no bearing to the problem because the question when 
a given particle arrives makes sense only when the parti- 
cle is in motion. We expect then that T is well defined if 
the contribution of the null subspace is removed. 

Technically this can be accomplished as follows [3Ï| : 
Let E and be the projections unto the closures of 
the subspaces Aí(p) (the subspace spanned by the zero 
momentum state) and ■Aí{p) ± (the subspace spanned by 
the non- vanishing momentum states), respectively. Now 
Aí(p) and A/'(p) ± are invariant under p; both subspaces 
then reduce p. Because p is self-adjoint, its restrictions 
on Aí(p) and Af(p) ± are both self-adjoint. The restric- 
tion pj^j-± has trivial null-space, so that its inverse, P^i, 
exists in r P ± TÍ. But this inverse exists only in V^Ti and 
not in Tí. This can be addressed by extending p7r± in 
the entire Tí. First we note that pT. 1 ^ is self-adjoint and 
bounded. It can be shown that pT,- 1 ^ is subnormal and it 
admits a unique minimal extension in the entire Tí. Its 
minimal extension is the bounded and self-adjoint oper- 
ator P^ 1 = E ± pJ^ ± E ± . This operator can be interpreted 
as the quantization of the classical observable for 
p =/= under the boundary condition imposed upon the 
momentum operator. 

Substituting P _1 for p _1 in the formal time of arrival 
operator, we get the operator quantization of T for peri- 
òdic boundary condition, 



1 



T = -- / i(qP- 1 + p- 1 q). 



(7) 



Note that both q and P _1 are bounded, everywhere de- 
fined, self-adjoint operators. Since T is symmetric under 
the exchange of q and P _1 , it is likewise bounded, every- 
where defined and self-adjoint. We shall refer to To as 
the periòdic confined quantum time of arrival operator. 
In position representation, To likewise assumes the form 
of a Fredholm integral operator, 



(TovXq) = / (q\T \q')v(q')dq'. 



-1 



whose kernel is 



q\ T oW) = 7-r( < ?+9')sg n (9-'?')- 



'2\ 



AifÜ 



(8) 



(9) 



The kernel (q\ Tq \q'} is likewise symmetric and square 
integrable. This means that the finite periòdic limit of 
© generates a self-adjoint integral operator, To, whose 
kernel is given by equation (JÜJ) . This operator is likewise 
compact — and thus discrete. And its eigenfunctions form 
a complete set of orthonormal system. 
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D. The Noncovariance of the CTOA-operators 

A time operator T is covariant if its eigenvectors, |r) , 
satisfy the property e - lHt / h | T ) — | T — t) , where \t — t) 
is also an eigenvector of T for any time t. Covariance 
of T implies that it has a completely continuous spec- 
trum taking vàlues in the cntirc real line. The time of 
arrival operator in L 2 (— oo, oo), for example, is covariant. 
Since the confined time of arrival operators posses a pure 
point spectrum, they are not covariant. Covariance has 
been a premium req uirement imposed upon time opera- 
tors [ïHlïlIïlIïlHElEllil· We wiU, however, 
demonstrate in the following sections that non-covariant 
time operators are physically meaningful. 



IV. THE CONJUGACY OF THE 
CTOA-OPERATORS WITH THEIR 
HAMILTONIANS 

A. Non-periodic case 

The commutator between T 7 and H 7 in the system 
Hilbert space is defined only if there exists a non-trivial 
intersection between the domains of the the composition 
operators T 7 H 7 and H 7 T 7 . It is not necessary that the 
commutator domain — the subspace T> com in which the 
operator (T 7 H 7 — H 7 T 7 ) is defined — coincides with the 
canonical domain — the subspace T> can in which the op- 
erator (T 7 H 7 — H 7 T 7 ) is proportional to the identity op- 
erator in T> C an- Generally we have the inclusion relation 
Dean Ç T> com for any pair of operators. In our case, we 
will find that T> can is a proper subspace of T> com . 



First let us find the domain of H 7 T 7 . 



Since T 7 is 



bounded, and thus defined in the entire Hilbert space, 
the domain of H 7 T 7 consists of all vectors <fi in the 
Hilbert such that T 7 is in the domain of the Hamil- 
tonian. Recali that the domain of H 7 consists of those 
that satisfy some boundary conditions (see below). Let 

4fism7 J_i 



-TT—z- %1 \ (l + q'Mq'W (10) 
4fism7 J q 

If tp(q) were to be in the domain of the Hamiltonian, 
first it must satisfy the boundary condition <p{—ï) = 
e~ 2l·1 <p{l) . Evaluating eauation- UlOJI at the boundaries 
yield 



<P{1) 



ne 11 í l 



Ah sin 7 



(/ + q')<t>(q')dq' 



<p(-l) = - 



//e 



-17 



Ah sin 7 



(-1 + q')0(q')dq' 



Imposing the first boundary condition on ip, gives us the 
equality J, (p(q')dq' = — J_, (j}(q')dq' ', which is only true 



if and only if both sides are equal to zero. Then <fi must 
satisfy the condition J , <f>(q')dq' = 0. That is the domain 
of H 7 T 7 is orthogonal to the subspace spanned by the 
zero-momentum-eigenfunction. 

Moreover, ip(q) = ÇT^4>)(q) must satisfy the second 
boundary condition <p'{—ï) = e~ 2n (p' (l) . Taking the first 
derivative of eauation- UlUfl gives 



V'(q) = - 



Ah sin 7 



' f%(q')dq' + e- 1 ^ f cj)(q')dq' 



(11) 



Because <j>{q) must satisfy J_ t (f){q')dq' = 0, the vàlues 
at the boundaries of the derivative simplifies to <p'{ï) — 
—i^jih~ l l(j)(l) and tp'(-l) — i^ih^l^—l)- Imposing the 
second boundary condition gives us 4>{—l) = —e~ 2l ' y (j)(l). 
Note that there infinitely many vectors satisfying these 
in the domain of the Hamiltonian, and vectors that lie 
outside the domain of H 7 . 

The domain of the operator H 7 T 7 then consists of all 
vectors <j>{q) in the Hilbert space satisfying the condi- 
tions f^fiiq^dq' = and 4>{—l) = — e~ 2%1 4>{l). Because 
of the first condition, the domain is orthogonal to the 
one dimensional subspace spanned by the zero momen- 
tum eigenvector. The operator H 7 T 7 is then not densely 
defined. For all vectors cj> in this domain, H 7 T 7 acts as 



(H 7 T 7 0)(q) = -ih(j){q) + -ihq4>'{q). 



(12) 



We arrive at this expression by further diffcrentiating 
Eauation- (|ll·|) and by multiplying the appropriate con- 
stants. 

On the other hand, the domain of T 7 H 7 consists of all 
vectors ip in the domain of H 7 such that Hjtp is in the 
domain of T 7 . But since T 7 is bounded, taking the entire 
Hilbert space as its domain, the vector H 7 </5 is automat- 
ically in the domain of T 7 . The domain of T 7 H 7 is then 
the entire domain of the Hamiltonian. In this case T 7 H 7 
is densely defined because the domain of the Hamiltonian 
is dense. In this domain, the operator T 7 H 7 acts as 

(T 7 H 7 0)( 9 ) = -i»^(g) + iW(g) + i«e i T0'(-O,(13) 

where two successive integration by parts have been made 
to arrive at this expression, and the boundary condition 
on the elements of the domain of H 7 has been imposed 
in the simplification. 

Now the commutator [T 7 , H 7 ] = T 7 H 7 — H 7 T 7 is de- 
fined only on the subspace of the Hilbert space which is 
the intersection of the domains of the operators H 7 T 7 



and T 7 H 7 . 



The vectors in the domain of H 7 T 7 sat- 



isfy the boundary condition 4>(—l) = — e 2n (j>(l); while 



those in the domain of T 7 H 7 satisfy 4>{—l) = e 2lJ (j){l) 



-l) = e- 2t ^ 

In order for these two boundary conditions to be sat- 
isfied simultaneously, we must have <j>(—l) — 4*{l) = 0. 
Then the commutator domain consists of all vectors, ó 
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in the domain of the Hamiltonian satisfying the condi- 
tions f_ t <t>{q')dq' = and <j>(l) = <p(-l) = 0. In this 
domain, the commutator of H 7 and T 7 is 

((H 7 T 7 - T 7 H T )0)(g) = ihcj>(q) + I«e*V(-0. ( 14 ) 

Clearly H 7 and T 7 are not canonically conjugate in the 
entire commutator domain. However, restricting the do- 
main to those whose first derivatives vanish at the bound- 
aries gives us a canonical domain. Thus in the subspace 
of the domain of the Hamiltonian consisting of all vectors 
4>(q) satisfying 

j 4>{q')dq' = 0, 4> {k) {±l) = for k = 0, 1, (15) 

the Hamiltonian and the confined time of arrival operator 
are canonically conjugate, 

((H 7 T 7 -T 7 H 7 )0)(g)= ï ft0(g), (16) 

and they are conjugate in a non-dense subspace, which 
is not the usual for canonical pairs. 

B. Periòdic case 

Following the same steps above, we find that Ho and 
To form a canonical pair in a non-dense subspace of the 
Hilbert space consisting of the vectors satisfying the con- 
ditions 

ip [k) (±1) = 0, J q k (p(q) dq = for k = 0, 1. (17) 
That is 

((H T - T HoM?) = *'M«) (18) 

for all <p in the canonical domain. Being orthogonal to 
the two-dimensional subspace whose elements are ip(q) = 
a + bq for complex a and ò, the canonical domain is not 
dense. As in the former case, the canonical domain is 
smaller than the commutator domain of Hq and Tq. 



for the pair (H 7 ,T 7 ), one may question whether H 7 and 
T 7 can be appropriately labeled as a quantum canonical 
pair. A detailed answer to this issue has already been 
given by one of us in Reference-|40j. with which we refer 
the reader to. 

It is sufficient to point out here that the canonical com- 
mutation relation (CCR) [Q, P](p = ihip possesses nu- 
merous non-unitary equivalent solutions in a separable 
Hilbert space, and one such solution is the pair (H 7 , T 7 ). 
The set of properties of a specific solution is conseqüent 
to a set of underlying fundamental properties of the sys- 
tem under consideration or to the bàsic definitions of the 
operators involved or to somc fundamental axioms of the 
theory or to some postulated properties of the physical 
universe, so that there is no preferred solution to the 
CCR. 



V. SYMMETRIES OF THE TOA-OPERATORS 
AND RELATIONS AMONG THEIR 
EIGENFUNCTIONS 

In this section we derive the symmetries of the con- 
fined time of arrival operators, and from these symme- 
tries we derive the bàsic properties of their eigenfunctions 
and eigenvalues. And from these symmetries we will in- 
fer the relationships among the eigenfunctions and eigen- 
values for different vàlues of the boundary parameter 7. 
Most important is the identification of these symmetries 
as anàlogues of the classical symmetries of the classical 
time of arrival. 

Central to our discussion are the behaviors of the time 
of arrival operators under parity, l·l, and under time re- 
versal, 6, operations. Both operators are bounded and 
act on all vectors of the Hilbert space with the follow- 
ing corresponding operations in coordinate representa- 
tion n</j(ç,í) = íp(—q,t) and Qcp(q,t) = tp*(q,—t), re- 
spectively. In momentum representation, the actions of 
the parity and the time reversal operator are Hip(k, t) = 
ip(—k,t) and Qip(k,t) = <p*(—k,—t), respectively. In the 
following discussions, we will only consider the vectors at 
t = so that reference to the parametric time t can be 
omitted. 



C. Quantum canonical pairs 

The above prescribed quantization of the Hamiltonian 
and the time of arrival for the spatially confined particle 
yields the correspondence 

{T,H} = 1 -> [H 7 ,T 7 ] c il 7 , (19) 

where l 7 is the identity in the closure of the canonical do- 
main T>J an . In most cases of quantum canonical pairs, the 
commutator and the canonical domains coincide and the 
canonical domain is dense. Because these domains do not 
coincide and that the canonical domain being non-dense 



A. Non-Periodic 7 7^ ■§ Case 

The symmetries of the non-periodic TOA-operators 
follow directly from the invariance of their kernel under 
the following operations, 

{q\J- t \q') = -(-q\T,\~qy (20) 

<ç|T 7 |g') = -<g|T_ 7 | g y (21) 
<ç|T 7 |g') = (- q \T_ 7 \- q >) (22) 

We will find below that the above properties of the kernel 
dictates the properties of the eigenfunctions and eigen- 
values of the T 7 's among themselves. 
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1. Symmetry-1 

Let us derive the symmetry arising from equation- l|20[l . 
Let tp be any vector in the domain of T 7 , which is the 
entire Hilbert space, then 

0»(ç)= J {q\T^\q')ip{q'W (23) 

Acting both sides of this equation by the parity operator 
and changing variables in the integration by q' — > —q', 
we arrive at 

(ITT»(ç) = j (-q\ T 7 \-q') (U<p)(q')dq' (24) 

Acting both sides of this equation by the time reversal 
operator yields, 

(enT 7 ^(ç) = J {-qlT^-q')* (QU V )(q')dq' (25) 

where the identity (QHp)(q) = tp*(—q) has been used. 
Applying equation- J5DJ) finally gives (TjQHip)(q) — 
— (<dHT 7 <p)(q). Since this relation holds in the entire 
Hilbert space, we get the following combined parity and 
time reversal symmetry of the T 7 , 

n- 1 e^ 1 T 7 en = -t 7 . (26) 

From equation l|26[l we can infcr the relationship 
among the eigenvalues and eigenfunctions of T 7 for a 
fixed 7. We know that T 7 is self-adjoint and com- 
pact for a given 7, and thus the eigenvalues are real 
and countable, in particular, they are either positive or 
negative. Let f a ,-y be an eigenfunction of T 7 with the 
corresponding eigenvalue r CT . 7 7^ 0, where a constitutes 
the collection of quantum numbers necessary in specify- 
ing the eigenfunctions of T 7 . Using equation we 
have T 7 <^ CTi7 = — n _1 _1 T 7 9ri(^ (Tj7 . Because ip a ^ is an 
eigenfunction of T 7 with the eigenvalue r CTi7 , we get the 
relationship — T CTi7 yv i7 — II T-yQn^,- 7 . And this 
implies the eigenvalue relation 

T 7 0n</Jo- j7 = — Ta^QILipa^. (27) 

Thus QIÍip an is an eigenfunction of T 7 with the eigen- 
value —T a ,i- Since t CT)7 is not zero, the eigenvalues of 
íp ar1 and 8II^ ffj7 have equal magnitudes but with oppo- 
site signs. We arrived at this conclusion from the reality 
of the eigenvalue. 

Thus we have identified one quantum number s which 
takes on either ±1, indicating the sign of the eigenvalue. 
Later on we will find that s is related with the direc- 
tion of propagation of the eigenfunctions. We indicate 
s by writing the eigenfunctions in the form ^ 7 j where 
the (+)-sign indicates that it corresponds to the positive 
eigenvalue, and the (— )-sign indicates that it corresponds 
to the negative eigenvalue. We can now write also their 



corresponding eigenvalues as T^„. In particular, we have 
the relationships 

^- 7 = 6LV+ 7 r- 7 = -r+ 7 (28) 

where a now constitutes the rest of quantum numbers 
less s. Thus for every 7 and <r, there corresponds two 
eigenfunctions <p^ n which are related according to equa- 
tion (mi . 

In position and momentum representations, equation 
(|28ü leads to the eigenfunction relationships ip~ 7 (q) = 

<PÍa(-q) and <Pa n ( k ) = <Pa n ( k ), wher e ft n {k) = 

í-l < ^fe 7 ' ) * (i) dq, in which k takes the discrete vàl- 
ues k — 0, ±1, ±2, .... And these lead to the probability 
density relations 

Wt n {q)\ 2 = |vç. 7 (-9)| 2 . kí. 7 ( fc )| 2 = |<^, 7 ( fc )| 2 ■ 

(29) 

Equations (|29[) mean that the position distributions cor- 
responding to tp~ and ^5+ are mirror images of each 
other, and the momentum distributions corresponding 
to the same eigenfunctions overlap. 

2. Symmetry-2 

Also using symmetry (121(1 one can show that 
(6T 7 i/?)(ç) = — ÇT^ 1 Qíp)(q) for all ip in the Hilbert 
space. This leads to the symmetry relation 

e _1 T_ 7 e = -T 7 . (30) 

Now let <y9^ 7 be the eigenfunctions of T 7 for a given a. 
Then, with our established notation above, we have the 
following relationship 

= e ^7' T Í-7 = ~ T ^7- ( 31 ) 

In position and momentum representations, equation 
(|31|l leads to the eigenfunction relations (f^ _ Aq) — 

<Pv* 7 (q) and ip^_ 7 (k) = And these lead to 

the probability density relations 

|<y 5 í- 7 (9)| 2 = |^. 7 (<?)| 2 > |<^,- 7 (fc)| 2 = Wa n {~ k )\ 2 ■ 

(32) 

That is, the position distributions corresponding to (p^ 

and <p~ 7 overlap; and the momentum distributions cor- 
responding to y>+ 7 and ip~„ are mirror images of each 
other. 

3. Symmetry-3 

Using symmetry (1221 it can be shown that 
(T 7 Híp)(q) — (HT 1 (p)(q) for all ip in the Hilbert space. 
This implies the symmetry relation 

n _1 T_ 7 n = T 7 . (33) 
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Given the eigenfunctions of T 7 , (/?^ 7 , and the symmetry 
relation (|33[) . it can be shown the eigenfunctions of T_ 7 
are also given by 

^_ 7 = n V ± 7 , r <T ± _ 7 =r± 7 (34) 

These give the eigenfunction relations ip a „(q) = 
ip^(-q) and ^_ 7 (fc) = (p* y (—k). From this follows 
the probability relation 

kí,- 7 (<?)| 2 = Wt n i-q)\ 2 > |<^,_ 7 (*o| = |<^ 7 (- fc )| 2 

(35) 

These mean that the position and momentum distribu- 
tions corresponding to ip<r,-i an< 4 fcr,-y with the same 
eigenvalues are mirror images of each others. 

Below the above symmetries and probability relation- 
ships will be identified with the different symmetries of 
the classical time of arrival. 

B. Non-Periodic 7 = § and Periòdic 7 = Cases 

For 7 = § and 7 = 0, we find similar behaviors. For 
both cases, the kernel (g|T 7 |g') is invariant under the 
following operations, 

(qlTjW) = -<g|T 7 |g'>* (36) 
(<z|T 7 |g'> = (-q\T^\-q') (37) 

These symmetries dictate the relationship among the 
eigenfunctions of the time of arrival operators for 7 — 
- 

u, 2 . 

Following the same method employed above, equations 
()36J) and (|37J) lead to the following symmetries of the time 
of arrival operators for 7 = 0, 

9" 1 T 7 e = -T 7 (38) 

n^ 1 T 7 n = T 7 (39) 

Likewise using the same arguments used above, equation 
(|38|l leads to a pair of eigenfunctions with equal magni- 
tudes of eigenvalues but with opposite signs, i.e. <^jr 
and t^ 7 . In particular equation (|38|l yields the following 
relationships between the eigenfunctions corresponding 
to the positive and negative eigenfunctions, 

¥Ç,7 = e PÍ, 7 > T *.~/ = - T t, r (4°) 

On the other, hand equation (|39|) implics that the eigen- 
functions <^ 7 are likewise eigenfunctions of the parity 
operator with even parity. Thus we have symmetric po- 
sition and momentum distributions, i.e. 

|< 7 (g)| 2 = |^ 7 (g)| 2 , |< 7 (fc)| 2 = K#)| 2 (41) 

ovcrlapping both in position and momentum represen- 
tations and symmetric about the origin. These imply 
the position and momentum operators have zero expec- 
tation vàlues for these two cases, which further imply 
that the eigenfunctions for these two cases correspond 
to the classically indeterminate case of vanishing initial 
position and momentum. 



C. The Classical and Quantum Symmetries 

Earlier we raised the question which of the many time 
of arrival operators to consider; in particular, which of 
the T 7 is the appropriate time of arrival operator. We 
assert that all should be taken into account. We note that 
the classical time of arrival operator satisfies the follow- 
ing symmetries: t(q,p) = — í(— q,p), t(q,p) = — t(q, — p), 
and t(q,p) = t(—q,—p). Comparing these relationships 
with equations (|29|) . (|32[) . and we find that there is 
a perfect correspondences between these sets. In partic- 
ular, we have the following correspondences 

t(q, P ) = -t(-q,p) <— > T 7 = -e- 1 n- 1 T 7 ne,(42) 

t(q,p) = -t(q,-p) «— > T 7 = -9- 1 T_ 7 9, (43) 
t(q,p)=t(-q,-p) «— » T 7 = n- 1 T_ 7 n. (44) 

It is then clear that all vàlues for 7 7^ should be ac- 
counted for in order to accommodate the entire symme- 
try of the classical time of arrival. Thus not a single value 
of 7 is a sufficient quantization of the confined classical 
time of arrival. Note that from the symmetry of the 
eigenfunctions for 7 = 0, w/2, the CTOA-operators T 
and Tjl correspond to the classically indeterminate case 
for q = and p = 0. 



VI. THE CTOA-EIGENVALUE PROBLEM 

In the previous section, we arrive at the symmetries of 
the confined time of arrival operators, and these symme- 
tries connect the different eigenfunctions and their cor- 
responding eigenvalues. But that is how far the symme- 
tries can give us. In this section, we solve the eigenvalue 
problem for the CTOA-operators, which is Fredholm in- 
tegral operator problem of the second type. Our method 
is to convert the integral equation into its (integro) dif- 
ferential form. We are going to exploit the symmetry 
properties that we have derived above. It will be suffi- 
cient for us to solve explicitly for the positive-eigenvalue- 
eigenfunctions, and from them derive their correspond- 
ing negative-eigenvalue-eigenfunctions via their symme- 
try relationship established in the previous section. 



A. Non-periodic Confined Time of Arrival 
Operators, 7 / 

Our problem now is to solve for the eigenvalue prob- 
lem T 7 <y9 r = Tip T , for the eigenfunction lç t and the cor- 
responding eigenvalue T,for positive r. The eigenvalue 
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equation is explicitly given by 



(q\T 7 \q')<p T (<l')dq' 



-/ 



f/e 



4fi,sin7 J _i 



(q + q')íç T {q')dq' 



Ah sin 7 



(q + q')<p r (q')dq', (45) 



Differentiating eauation- (|45(l twice using Leibniz rule of 
difïerentiating an integral, yields the following differential 
equation for the eigenfunction, 

d 2 (p T (q) jJiiqdip T {q) 3/ii 



dq 2 



rh dq 2rh 



+ ^r(?)=0. (46) 



The eigenfunctions are distinguished among the solutions 
of this differential equation by extracting and imposing 
the boundary condition from equation- l|45|l itself. Evalu- 
ating equation- l|45|l at the boundaries yields the following 
integro-boundary conditions on the eigenfunctions 



fi e 



«7 



4r?ïsin7 J_ 



\i e 



-17 



Arh sin 7 



(i + q ')<p T (q')dq', (47) 



(-l + q')<p T (q')dq'. (48) 



These are non-standard boundary conditions. Neverthe- 
less they are sufficient to determine the eigenfunctions 
and eigenvalues. 



1. The general case 

Now we solve differential equation (|46|l for a given 7 
subject to the conditions and (gHJ). Equation- (|4*6*)) 
has a deffnite parity. If (p(q) is a solution, then Lp{—q) 
is also solution. This can be seen by making the sub- 
stitution q — > — q in the differential equation. It is then 
sufficient for us to find odd and even solutions, and from 
them built the general solution by linear superposition. 
By power series method, we find the following odd and 
even solutions 



<Pe(q) 



<Po{q) = qe 4Ta 



The eigenfunctions are then of the form, 
íf T = A ip e (q) + AiLp (q) 



(49) 



(50) 



(51) 



where Aq and A\ are constants yet to be determined 
from the boundary conditions. Substituting equation- 
()51J) back into both sides of equations Ij47(l and (|48|l . and 
performing the indicated integrations, we obtain, after 
some simplification, the following system of equations for 
the unknown coefficients Aq and A\ , 



J 



A „M íj_,_ (£) + _L A (£)) + All fj_, (Q + _!_,, (£))_ 

V t% \ 4 \ÍTh J tan 7 4 \4ThJJ \ 4 \4ThJ tan7 4 \4t7ï / / 

V rh \ -\4tHJ tan7 ^ \AtH J J \ ^ \AtH J t&nj * \4tH J J 
This system of equations can be written in matrix form: 

/^V^(M^) +C0t ^(^)) r(!)(^(^) + cot7J f (^)) W A \ 



In order for a non-trivial solution to exist, the determi- 
nant of the matrix of the coefficients of Aq and IA\ must 
vanish. The vanishing determinant leads to the condition 

J-éx) J-í(x) - cot 2 7Ja(x) Ji(x) = 0. (52) 

where x — fiP/Arh. By functional analytic arguments, 
the roots of this equation must exist and real and count- 
ably many. 



If we order the roots of eouation- (|52|) such that n = 1 
corresponds to the first positive root ri, n = 2 to the 
second root ri-, and n — m to the m-th root r m , then 
we find that the remaining quantum number completely 
specifying the eigenfunctions of the confined non-periodic 
time of arrival operators consists of the positive integers 
n ordering the roots of the eauation- (|52|l . For a given 
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root r n , we find that the constant A\ has the form 



eigenvalue eigenfunction (g). 



M = y^( 3 -^\-^fA M (53) 

Substituting Ai back into equation- JSTJ gives the positive 
I 



And given tp+ (q) we can determino ip n (q) from the 

symmetry relation i.e. (Pn-y(l) = ®ïl<Pn -yil)- After 
pcrforming some simplifications, the eigenfunctions are 
given by 



+ 



(<?) 



2g 



i 6 " r " P 



2 \ 4 



J_ i(r„ ) - cot 7 J|(r„) J 

M_|(r„) - cot 7Ji(r n )) 



(54) 



where A na is the normalization constant. The corre- 



sponding eigenvalues are 



n,7 



±- 



'4fi r„ 



(55) 



We shall call those that do not vanish elsewhere in the 
interval [—1, l] as nonnodal-eigenfunctions; otherwise, as 
nodal eigenfunctions. The non-nodal (nodal) eigenfunc- 
tions correspond to those with even (odd) quantum num- 
ber n. 



2. The special anti-symmetric case, 7^0 

The eigenfunctions for T» can be derived directly 
from cauation- 15411 . For 7 = tt/2, eauation- (|52|l re- 

duces to J_dx) J_iix) = 0. Then either J_3(x) = or 

4 S ' 4 V ' 4 V ' 

J-x(x) — 0. For the later the second term of the eigen- 

4 

function given by equation- vanishes; for the former, 
on the other hand, the first term vanishes. In this case, 
the eigenfunctions bifurcate into odd and even eigenfunc- 
tions. The even and non-nodal eigenfunctions are 



J_ 



*P 
(56) 

with the eigenvalues given by rf* e = ±(^il 2 /4h)r s 1 , and 
the r s 's are the positive roots of J_a(x) = 0, with s 
1.2 And the odd and nodal eigenfunctions are 



<*,o(?)=4Uïe 



[ J -^4) TlJ f( r 4 

(57) 

with the eigenvalues given by t^ 2Lo = ±(pl 2 /4:h)r u 1 , 
and the r u 's are the positive roots of J_i(x) — 0, with 
u = 1, 2, . . . . A e s and A° u are normalization constants. 



B. Periòdic CTOA-Operator, 7 = 

We turn to the eigenvalue problem T </5 r = Tip T , for 
positive r. Explicitly, the eigenvalue equation is given by 

/q pi 
ip T (q + q')dq' - — J tp T {q + q')dq' 



l 



ip T (q -q' )dq' 



(58) 



Differentiating this equation twice and after performing 
some simplifications yield the integro-differential equa- 
tion 



d 2 ip T [iiq dip T 



dq 2 



T n h dq 2rh 



fií 
2TÜ 



-/ 



<p T (q')dq' (59) 



And the eigenfunctions are subject to the integro- 
boundary conditions 



ip T (l) 



Airh 



1 f 

(l + q')(p T (q')dq'—r / (l 2 -q' 2 )(p T (q')dq' 
l 1 J-i 



Airh 



(l-q'^Aq'W-j (l 2 -q' 2 )<pAq')dq' 



where we have arrived at these by evaluating equation- 
(|ïïü|) at the boundaries. 

For us to have an idea of the possible solutions of the 
integro-differential equation satisfying the above bound- 
ary conditions, we subtract and add the boundary con- 
ditions to yield 



tp T (l) - <f T (-l) 



tp T (l) + <Pr{—l) 



í2tH 



q'<p T (q')dq'. (60) 



-/ 



I' 



%21tK 



q' 2 Lp T (q')dq (61) 
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In this form of the boundary conditions, it is evident that 
solutions are odd and even functions of position. The 
even (odd) solution, for example, trivially satisfies the 
first (second), whilc it must satisfy the nontrivial second 
(first) condition. This is in fact predicted by our sym- 
metry analysis, where we found that the eigenfunctions 
eigenfunctions of the parity operator as well. 



1. The Odd Eigenfunctions 

For odd eigenfunctions the right hand side of equation- 
(|59|l vanishes and we are left with the same differential 
equation that we have for the non-periodic case. We 
know already that the odd solutions are of the form, 
which are given by eouation- (|50[l . For this case the 
integro-boundary condition l|60[) reduces to 



vo, (0 = -^y vo(q')q'dq'. 



(62) 



Substituting the odd solution given by eauation- ll5()l) into 
the right hand side of this boundary condition gives 
J_i(x) = with x = fil 2 /4Th. 

Now let r„ be the nth positive root of J_ i (x) . Then 
the positive eigenvalue corresponding to this root is 
T n = il· t ^ 2 /^) r n 1 - Substituting this eigenvalue back 
to eauation- (|ïïül) gives the corresponding odd positive- 
eigenvalue-eigenfunction <p„ )0iO (<7). Given ¥>^, 0o (?) 
we likewise have the negative-eigenvalue-eigenfunction 
fnfitoiq) usm g symmetry relation (|JjJ. i.e. <p n o (q) = 
®(p^ 0o (q). Explicitly, the odd eigenfunctions are given 
l'.v 



t fÍ,oAq)= A °nflq eTtru ^ T [ 



2\ 4 



with the corresponding eigcnvalucs 

fil 2 



rt = ±- 



(63) 
(64) 



~4r„/i' 

for n = 1,2,.... A° n is the normalization constant. 

2. The Even Eigenfunctions 

For even eigenfunctions, the right hand side of 
equation- (|59|l does not vanish. Since the right hand side 
involves only an integral of the unknown eigenfunction, 
the integral can be equated to a constant A, which is 
to be determined. This reduces the integro-differential 
equation into the form 



d 2 (po(q) + fnq dtpojq) + 3fii _ \ii 

dq 2 T n H dq 2rh 2tM 



A. 



(65) 



with the solutions subject to the integro-boundary con- 
dition 



li 



f<p T {q')dq'. 



(66) 



-/ 



Equation- (|65() is a linear, second - order non homoge- 
neous differential equation. To solve it, it is sufficient to 
find the general solution to its associated homogeneous 
differential equation, which happens to be just equation- 
14*51 If (^0(9) is the general solution to the associated 
differential equation, the general solution to l|65|l is given 
by 



^ 0) (q)^^(q) + f l A 



1 



(67) 



The even eigenfunctions then are found by obtaining the 
even solutions of the associated differential equation. 

But we know already the solution to cquation- (|46|l . 
Substituting the even solution (|4*9"|> back in equation |J(77|) 
and evaluating the solution at the boundaries yield the 
constant 



A = 12Zexp 



liil 2 



Ji 



\AtH 



(68) 



Substituting igni) and into JHII), we obtain the follow- 
ing analytic form of the solution to the integro-diffcrcntial 
equation, 



<A), e (<7) 



rh 



■ 7 ( m 

4 \4r/ï 



2 \ 1 



-4e 



~ 4i-h 



Ja 



Arh) 



(69) 



The eigenvalues are found by imposing the integro- 
boundary condition on these solutions. Substituting 
eouation- (|69|l back into (|66|l leads to the equality 

2 1 
J_a(x) + -Js(x) + -Ji(x) = 0. (70) 

4 w 34 w x 

where x = fil 2 /4rH. 

The eigenvalue problem is then reduced to finding the 
roots of eQuatioml70l Let r s be the sth positive root of 
this equation, the roots being ordered according to in- 
creasing magnitude. The s-th positive eigenvalue is then 
given by rf — . Substituting this back into equation- 
(I69|l gives the corresponding eigenfunction. From these 
positive eigenvalue eigenvalue eigenfunctions we derive 
the negative eigenvalue eigenfunctions using symmetry- 
(14011 . The eigenfunctions are now explicitly given by 



•/_4l y^TijJ^T, 



4^ ;0 e^Ji(r s ) 
(4r s )ï 



(71) 
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where A St o is the normalization constant. The corre- 
sponding eigenvalues are 



= ± 



fil 2 
4r.fi 



s = 0,±1,±2· 



(72) 



VII. DYNAMICS 

The question now arises as to how we should intèrpret 
the eigenfunctions and eigenvalues of the confined time of 
arrival operators. Standard quantum mechanics postu- 
lates that the eigenvalues of an observable are the results 
of measurements when that observable is subject to mea- 
surement. In this section, we will demonstrate numeri- 
cally that the spectral properties of these operators ac- 
quire interpretation independent from the measurement 
postulate. Here we will find that the spectral properties 
are instead tied with the dynamics of the system. 

From the symmetries of the CTOA-operators, it is ev- 
ident that the negative eigenvalue-eigenfunctions have 
exactly the same dynamics as those of the positivc 
eigenvalue-eigenfunctions in the time reversed direction. 
It is then sufficient for us to consider in detail the dynam- 
ical behaviors of the positive eigenvalue eigenfunctions. 
We will classify the evolution according to whether the 
eigenfunction concerned is an eigcnfunction of parity op- 
erator or not. The parity eigenfunctions are those of the 
7 = 0, 7r/2 cases, and non-parity eigenfunctions other- 
wise. 



A. The quantum equation of motion 

Our analysis is based on the numerical integration of 
the evolution law ip(t) — e _íH '' t/ '' í y>(0) using spectral de- 
composition method, which, in position representation, 
is explicitly given by 



■iE k t/hil 



(q) 



e ikn J 



k— — ' — 1 



(73) 



where b k {t) = e - iEkt / h J l _ t 4>l*(q)<p(q,0)dq. Evidently 
the time evolution is a reconstruction by Fourier series. 

Central to our analysis for the interpretation of the 
eigenfunctions and eigenvalues are the expectation value 
and the variance of the position operator as a function 
of time with respect to the eigenfunctions of the confined 
time of arrival operators. That is the quantities 



q\(p nil {q,t)\ dq 



-/ 



(75) 

where the <p n ,y (g,í)'s are the evolved eigenfunctions of 
the confined time of arrival operators. 



B. Numerical evolution and the Gibbs phenomenon 

Since there is no closed form for the Fourier coefficients 
of the eigenfunctions, we resort to numerical evaluation 
of the coefficients and hence a numerical evaluation of the 
sum involved in the evaluation of eauation- H73(l . The nu- 
merical implementation of equation- JJ^Í is then a special 
case of the general truncated Fourier series 



TV 



Ín(x)= J2 fi 



k e 



ik-nx 



(76) 



k=-N 



where the fk's are the Fourier coefficients given by, in the 
rescaled interval [— 1, 1], fu = / , e~ tkvx f(x)dx. 

It is well known that when f(x) is analytic and peri- 
òdic the Fourier series converges exponentially fast, i.e. 
max \f(x) — fjsr(x)\ < e~ aN for — 1 < x < 1 for some 
a > 0. In such cases the truncated Fourier sum is 
an accurate approximation of f(x) for sufficiently large 
N. However, when f(x) is non-periodic and/or discon- 
tinuous functions the Fourier sum converges slowly in 
the interval (-1,1), i.e. \f(x ) - f N (x )\ » O^" 1 ) 
for — 1 < Xq < 1. And there is an overshoot at the 
boundary that does not diminish with increasing N, i.e. 
max \f(x) — //v(x)| for — 1 < x < 1 does not tend to zero 
as N increases indefinitely. This is the well-known Gibbs 
phenomenon, which undermines obtaining accurate point 
vàlues from the knowledge of Fourier coefficients for non- 
periodic functions |42|. 

The presence of Gibbs phenomenon undermines our in- 
tention to understand the dynamics of the eigenfunctions. 
This can be seen from eauations- (|74f) and-J7SJ where we 
need accurate point vàlues of the evolved eigenfunctions 
in order to get an accurate value of the required inte- 
grals. In order to have an accurate picture of the evolu- 
tion of the expectation value and variance of the position 
operator, one must have first an accurate numerical ap- 
proximation to the evolved eigenfunctions. Thus in the 
following we limit our numerical evolution to particular 
set of eigenfunctions and vàlues of 7 in order to avoid 
Gibb's phenomenon. 



Results 



1. Parity eigenfunctions 



(q)n,.W 



q IWyfà, í)| dq 



(74) 



The eigenfunctions for the symmetric and anti- 
symmetric confined time of arrival operators are parity 
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FIG. 1: The variance as a function of time for the first ten 
symmetric CTOA-even-eigenfunctions for 7 = 0, (a); and an- 
tisymmetric CTOA-odd-eigenfunctions for 7 = 7r/2, (fe). The 
corresponding minimum variance for these eigenfunctions are 
shown in(c). I = fj, = h = 1. All units, with all in succeding 
figures, are in atòmic units. 

eigenfunctions, being even and odd functions of posi- 
tion. For the 7 = case we can only evolve directly 
the even eigenfunctions via Fourier series without the 
effect of Gibbs phenomenon. For the 7 = tt/2 case, 
we can evolve odd eigenfunctions without the effect of 
Gibbs phenomenon. That that is the case for the anti- 
periodic case can be seen from the following argument. 
The Fourier coefficicnts are given by 

J_4*(l)Vn,%(q)dq ~ fe-^T (e^^n.f («)) dq. 

Evidently the sum in equation- J73J) is the Fourier sum for 
the analytic function (e -I7r ~ <p n ,%(q)) ■ For even ip n ^(q) 
the function (e~ l7T Tip n ^ (q)) is odd, thus nonperiodic; 
but for odd eigenfunction, it is even, so that it is periòdic. 
Here we demonstrate numerically that these eigenfunc- 



tions unitarily arrive, i.e. the eigenfunctins evolve ac- 
cording to Schrodinger's equation, at the origin at their 
respective eigenvalues. These unitary arrival is measured 
by the variance of the position operator with respect 
to these eigenfunctions. Figure-Q] shows the variance as 
a function of time for the first ten largest eigenvalue- 
eigenfunctions for the even and odd periòdic and anti- 
periodic CTOA-eigenfunctions, respectively. Starting 
from t = the variance decreases, reaches a minimum 
then increases again, followed by a decaying oscillation. 
The figure clearly dcmonstrates that the variance de- 
creases with n, so that the eigenfunctions are arbitrarily 
localized at the origin for arbitrarily large n. Tables- I 
and-II shows the minimum variance for a given eigenfunc- 
tion and the interpolated variance at the eigenfunction's 
corresponding eigenvalue. For the parameters indicated 
in the caption, the minimum variance and the interpo- 
lated variance at the eigenvalue agree at least to three 
significant figures. Our results shows that the eigenfunc- 
tions evolved such that the variance is minimum at their 
eigenvalues. 

Figure-|2| shows the general features of the evolving 
probability density for even and odd parity eigenfunc- 
tions. The even parity eigenfunctions evolve such that 
their corresponding probability densities obtain their 
minimum widths, and their peaks being maximum at 
the origin at their respective eigenvalues. The proba- 
bility densities for the odd eigenfunctions, on the other 
hand, evolve with two peaks approaching the origin, the 
value of the probability density being zero at the origin 
for all times. The time of closest approach to the origin 
of the two peaks oceur at the eigenvalue. We shall refer 
to the former as non-nodal and the latter nodal. We will 
find below that the dynamical behaviors of these eigen- 
functions are also shared by the eigenfunctions for the 
non-periodic CTO-operator eigenfunctions. 



2. Non-parity eigenfunctions 

The eigenfunctions for 7 =/= are non-parity eigenfunc- 
tions. These eigenfunctions are nonperiodic, hence sub- 
ject to Gibbs phenomenon. However, we can choose 7 
such that Gibbs phenomenon can be "reduced" . Notice 
that the the eigenfunctions are linear superpositions of 
odd and even functions of q. The idea is to choose 7 such 
that the even part dominates the odd part. If this can be 
done, the eigenfunctions are approximately periòdic, and 
hope that Gibbs phenomenon is not that "large" . Fortu- 
nately, this can be done. Recali that the eigenvalues are 
determined from the roots of the characteristic equation 
J- 3 (x)J_i (x) — cot 2 7J3 (x) Ji (x) = 0. For sufficiently 
small 7, the second term of the characteristic equation 
dominates and the first term and the eigenvalue condi- 
tion reduces to Js (x)Ji (x) ~ 0. The eigenvalues are now 

4 W 4 W 

approximated by the roots of Jz (x) and Ji {x) . For the 
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Periòdic Even CTOA-eigerifunctions 



n-th 
even 
c i ge ii- 
function 


Eigenvalue 
(xlO J 


Computed 
minimum 
variance in the 
interval 

/ v i n— 3\ 

(xiu ) 


Interpolated 
variance at the 
eigenvalue 


i 


111.43823 


122.4 


122.6 


2 


50.17966 


63.72 


63.76 


3 


31.58532 


39.72 


39.72 


4 


22.84895 


27.63 


27.63 


5 


17.84116 


20.76 


20.76 


6 


14.61367 


16.46 


16.46 


7 


12.36664 


13.57 


13.57 


8 


10.71462 


11.50 


11.50 


9 


9.44995 


9.970 


9.971 


10 


8.45118 


8.787 


8.787 



TABLE I: Comparison of the calculated minimum variance 
with that of the linearly interpolated variance at the eigen- 
value for the corresponding eigenfunction for the first ten even 
functions of the periòdic CTOA-operator. The variances are 
calculated with 401 Fourier coefncients at 0.0001 time steps. 



Anti-Periodic Odd CTOA-eigenfunctions 



n-th 


Eigenvalue 


Computed 


Interpolated 


odd 


(xlO" 3 ) 


minimum 


variance at the 


eigen- 




variance in the 


eigenvalue 


function 




interval 








(xlO -3 ) 


(xlO -3 ) 


1 


124.60751 


161.0 


161.3 


2 


48.79893 


52.46 


52.47 


3 


30.27385 


30.61 


30.61 


4 


21.93662 


21.50 


21.50 


5 


17.19832 


16.53 


16.53 


6 


14.14287 


13.42 


13.43 


7 


12.00910 


11.29 


11.29 


8 


10.43469 


9.741 


9.746 


9 


9.22520 


8.558 


8.558 


10 


8.26694 


7.631 


7.633 



TABLE II: Comparison of the calculated minimum variance 
with that of the linearly interpolated variance at the eigen- 
value for the corresponding eigenfunction for the first ten odd 
functions of the anti-periodic CTOA-operator. The variances 
are calculated with 401 Fourier coefncients at 0.0001 time 
steps, except for the 9-th and 10-th eigenfunctions where the 
time step is 0.00005. 

eigenfunctions corresponding to the roots of the 
even term dominates the odd term, and the eigenfunc- 
tions are approximately periòdic at the boundaries. On 
the other hand, for the eigenfunctions corresponding to 
the roots of J|(x), the eigenfunctions are non-periodic, 
and we expect that Gibbs phenomenon has considerable 
effect on the numerical sum. Figure-|3] gives a graphi- 
cal comparison of the The functions Ji(x), J^x), and 
J- 1 (%)J- i [%) — cot 2 7J3 (x) Ji (x) for 7 — 0.01. It shows 
that for the given 7 the roots of Ji (x) and J| (x) approx- 
imate the roots of J_ 3 (x) J_ 1 (x) — cot 2 7J3 (x) Ji (x). 




FIG. 2: (Color online) The n = 10 (a) evolved probability 
density for 7 = 0, and n — 10 (b) evolved probability density 
for 7 = 7r/2, with h = l = m = 1. Both symmetrically 
"collapse" at the origin at their respective eigenvalues. 
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FIG. 3: The functions Ji(x) (sòlid line), J3 (x) (dashed line), 

4 4 

and J_ 3 (x) J_i(x) — cot 2 (0.01) J3 (x) Ji (x) (the almost ver- 

44 44 

tical lines). 



Figure-Q]a shows the behavior of the expectation value 
of the position operator for each eigenfunction for 7 = 
0.01. We find that the eigenfunctions for even n approach 
the origin from the positive side axis; whilc the eigenfunc- 
tions for odd n, from the negative axis. Figure-^f) shows 
the variance as a function of time. Evidently the vari- 
ance has the same behavior as those of the parity eigen- 
functions for 7 = 7r/2,0; and that the minimum vari- 
ance decreases with n so that the eigenfunctions become 
increasing localized at the origin for increasing n. The 
eigenfunctions are also either nodal or nonnodal with the 
same dynamical behaviors as those of the parity eigen- 
functions. For the given 7, the nodals are those with n 
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0.01 5 





FIG. 4: The top-figure shows the expectation value of the 
position operator for the first ten (10) eigenfunctions for 7 = 
0.01. The sòlid line corresponds to n equals even; dashed line, 
n equals odd. The bottom-figure shows the variances, c 2 (í), 
of the same eigenfunctions as a function of time. I = /1 = h — 
1. 

odd; and non-nodals, with n even. Figurc-|S] shows the 
general features of the probability density as a function 
of time and space for nodal and non-nodal eigenfunc- 
tions. Generally the zero of the nodal eigenfunctions do 
not occur at the origin as those of the nodal parity eigen- 
functions; however, the zero approaches the origin in time 
and coincides with the origin at the eigenvalue. 

Table-III summarizes the minimum variance, together 
with the linearly interpolated position expectation value 
and variance at the eigenvalue. It shows that at the 
eigenvalue, the expectation value is qualified zero for all 
cases. On other hand, the minimum variance and the 
interpolated variance at the eigenvalue match at least to 
three significant figures for n-even and at most two sig- 
nificant figures for n-odd. The former correspond to the 
case when the eigenfunction is approximately periòdic at 
the boundary as discussed above; the latter, when the 
eigenfunction is non-periodic and Gibbs phenomenon is 
prevalent. The blank entries for the largest eigenvalue 
indicates the failure of our method to converge. This is 
due to the large oscillations in the exponentials for long 
times of evolution. 



D. The interpretation of the CTOA spectral 
properties 

The foregoing numerical results strongly endorse the 
minimal interpretation that a CTOA eigenfunction is 
a state that evolves to unitarily arrive at the origin at 
its eigenvalue — that is, a state that evolves according to 



FIG. 5: (Color online) The n = 20 (a) and n = 21 (b) evolved 
probability densities for 7 = 0.01, with h = l = m = 1. Both 
unitarily arrive at the origin at their respective eigenvalues, 
0.0081 and 0.0079. 



Non-Periodic CTOA-eigenfunctions 7=0.01 



n-th 


Eigenvalue 


Average 


Computed 


Interpolated 


Eigen- 


(xl0~ 2 ) 


Position 


minimum 


Variance 


function 




at the 


variance 


at the 






Eigenvalue 




Eigenvalue 






(xlO" 4 ) 


(xlO" 2 ) 


(xlO" 2 ) 


1 


2887 








2 


8.990 


-2.40 


9.629 


9.665 


3 


7.161 


1.19 


6.367 


7.208 


4 


4.233 


-0.49 


4.191 


4.192 


5 


3.758 


0.35 


3.361 


3.586 


6 


2.765 


-0.19 


2.653 


2.653 


7 


2.551 


0.16 


2.264 


2.366 


8 


2.052 


-0.11 


1.935 


1.936 


9 


1.931 


0.09 


1.697 


1.758 


10 


1.632 


-0.07 


1.522 


1.522 



TABLE III: Comparison of the calculated minimum variance 
with that of the linearly interpolated variance at the eigen- 
value for the corresponding eigenfunction for the first ten 
eigenfunctions functions of the periòdic CTOA-operator. The 
variances are calculated with 601 Fourier coefncients at 0.0001 
time steps. 



Schrodinger's equation such that the event of the position 
expectation value assumes the value zero, and the event 
of the position variance or uncertainty being minimum 
occur at the same instant of time equal to the correspond- 
ing eigenvalue. This provides the justification for the 
identification of the CTOA-operators as time of arrival 
operators. While our conclusion is based on a limitcd 
range of numerical results, these results are nevertheless 
compelling that we can take the minimal interpretation 
as an exact statement of the dynamical properties of the 
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eigenfunctions of the confined time of arrival operators, 
until otherwise is proven. 

But what is the importance of this realization? Recali 
that one of the surrounding issues in the past against 
the lcgitimacy of the quantum time of arrival problem 
within Standard quantum mechanics is the absence of 
phase space trajectory for a quantum particle, so that the 
question of quantum time of arrival is ill-defined. Our 
results here demonstrate that the QTOA-problem can 
be rephrased within the Hilbert space and the problem 
translates to finding states that unitarily arrive at a given 
point at a definite time. The QTOA-problem phrased 
in this way is well-defined because quantum states have 
well-defined trajectòries in the Hilbert space according to 
the Schrodinger equation. 

But does this not trivialize the QTOA-problem be- 
cause one can always construct by hand states that will 
evolve to have arbitrarily sharp width at the origin by 
exploiting the unitarity of quantum dynamics [43t l44j |? 
Constructing a single state with such properties is trivial, 
but the problem becomes non-trivial when we require the 
solution to comprise a complete set of such states, with 
the operator that can be constructed from this set being 
canonically conjugate with the Hamiltonian, and has an 
unambiguous classical limit which is the classical time of 
arrival 47]. These requirements maybe necessary if we 
were to preserve the quantum-classical correspondence. 
The completeness requirement is inevitable if we require 
a theory of quantum first time of arrivals that is reflec- 
tive of the fact that the entire phase space is accessible 
to a quantum particle via quantum tunneling. Clarifi- 
cation of these issues will have to wait for the general 
theory of confined quantum time of arrivals for arbitrary 
potentials, and the investigation of the limit as l goes to 
infinity llaliíi. 



VIII. CONCLUSION 

In this paper we have given full account of the con- 
fined quantum time of arrival operators. While our re- 
sults have answered some qüestions, they have raised 
some more qüestions and left others unanswcrcd. Now 
we know self-adjointness and conjugacy of a time oper- 
ator with a semi-bounded Hamiltonian can be achieved 
simultaneously; and a time operator need not be covari- 
ant and thus can be compact, with the eigenfunctions 
and eigenvalues tied with the dynamics of the system, 
acquiring interpretation independent from the quantum 
measurement postulate. 

However, these very realizations raise several funda- 
mental qüestions. The dynamical interpretation of the 
eigenvalues of the CTOA-operators does not fit well with 
the fundamental quantum measurement postulate for ob- 
servables represented by self-adjoint operators with dis- 
crete spectrum. Does this imply a reconsideration of the 
fundamental quantum measurement postulate? In par- 
ticular, does this call for a classification of observables 



according to their relationship or non relationship with 
the system dynamics? What sort of modification to the 
fundamental postulates to accommodate such classifica- 
tion if indeed necessary? But then this leaves us with the 
question of the relationship between the spectral proper- 
ties of a discrete time operator with the actual clicks of a 
detector. Should these qüestions prové imperative their 
answers will have non-trivial repercussions at the foun- 
dational level. 



Appendix A: Derivation of the kernels 

The kernel (q\ T 7 \q') for the non-periodic CTOA- 
operator is derived as follows. Using the property q | q) = 
q\q) , we have 



2 

JL 
Ah 



q + q')(q\p'. 



-1 l„'\ 



(q + q')e^-^ £ 



e i 



<q-q') 



7 + kn 



{V) 



where the second line follows from introducing a resolu- 
tion of unity provided by the complete eigenvectors of 
the momentum in the factor (q\ p^ 1 \q'). Now the sum 
can be rewritten in the form 



E 

k= — oo 



,-r i <9-*') [ 1 ^ cos [!f(q - q')] 



7 + kn 



k=l 



7 2 — fc 2 7T 2 



k=l 



7 2 — k 2 TT 2 



The two infinite series can be straightforwardly evaluated 
by contour integration. It yields the explícit forms 



E 

k=l 
oo 

E 

k=i 



COa[!f(q-q')] cos 
7 2 — fc 2 7T 2 

k s in[!f(q-q')]_ sin 



7 1 



\q-q' 



27 sin 7 



7 2 — fc 2 7T 2 



2-7T sin 7 



1 

sgn(ç - q'). 



Substituting these back into Equation- 1|78[) and after 
some simplification, we have 



E 

k— — oc 



7 + kir 



cos 



7 1 



\q~q'\ 
1 



sin 7 
7(l- 



\q-q'\ 
1 



sin 7 



■sgn(ç-g')- 



This can still be simplified, 



. (»-»') 



E 

k— — 00 



7 + k,7T sin 7 



cos 7, 



,(«'-») 



q > q 
q = q' 

q<q' 



(79) 
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Finally substituting this back into eauation- (|77(l give us 
the kernel- © , as long as we define H(0) = 1/2. 

The kernel (q\ Tq \ q') for the periòdic CTOA-operator 
can be derived similarly. Again using the property 

q I <?) =q\q) gíves us 



<?|Tok') 



= -^( q + q')(q\P- 1 W) 



JL 

Ah 

2nlh 



oc 1 



k=— oo 



oo ^ 



fc=l 



k-K , 

T (í-9 



Using the identity X^^Li n 1 sinir = 7rsgn(a;) — x for all 
— 2tt < x < 2tt, we finally arrive at the kernel-®. 



Appendix B: 
Numerical Solution to the CTOA Eigenvalue 
Problem 

An independent numerical solution can be obtained 
for the CTOA-eigenvalue problem. Here we describe the 
Nystrom method of solving the Fredholm integral eigen- 
value problem. Generally the integral operator eigen- 
value problem is of the form 



K(q,q')tp(q')dq' = \p(q), 



(80) 



where (p is an eigenfunction and A the corresponding 
eigenvalue, with a < b. The Nystrom method is based on 
some choice of an integration quadrature, f^?p(q)dq = 

YljLx w j' t l J (lj)í where the Wj's are the weights of the 
quadrature rule and the N points gy 's are the abscissa. 
Using the quadrature rule, eauation- ljSüf) reduces 
q^Wjifiqj) = \<p(q), and evaluating 
this cxpression at the absiccas further reduces to 



merically the Fredholdm eigenvalue problem then reduces 
to the matrix eigenvalue problem 



1,2, 



Nu- 



K ■ (p = Xíp 



(81) 



[—1, l] to the interval [—1, 1]. The N abscissas are the N 
roots of the Legendre polynomial 



P N (x) 



fc=0 



2N-2k 

N 



J-2k 



in the interval [—1, 1]. The weights are given by 



Wj 



j\P n (x) 2 dx 



(82) 



(83) 



where P' N (xj) is the derivative of the Legendre polyno- 
mial at its zero Xj. 

Table-IV shows the comparison of the exact and the 
numerical results for the eigenvalues of the CTOA- 
operators for the indicated vàlues of 7. The eigenvalues 
have been computed using Nystrom method and employ- 
ing Gauss-Legendre quadrature with 10000 integration 
points. There is an excellent agreement between the ex- 
act and the numerical vàlues for the CTOA-operators 
with definite parity, i.e. for 7 = 0, and a fair agree- 
ment for other vàlues of 7. 
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7 = 


7=£ 


7=z 


exact 


numèric 


exact 


numèric 


exact 


numèric 


1 


0.124608 


0.124608 


0.73965 


0.73953 


0.37834 


0.37829 


2 


0.111438 


0.111438 


0.09459 


0.09468 


0.10510 


0.10515 


3 


0.050180 


0.050180 


0.06887 


0.06878 


0.06415 


0.06411 


4 


0.048799 


0.048799 


0.04332 


0.04341 


0.04539 


0.04543 


5 


0.031585 


0.031585 


0.03683 


0.03674 


0.03545 


0.03541 


6 


0.030274 


0.030274 


0.02806 


0.02815 


0.02892 


0.02896 


7 


0.022849 


0.022849 


0.02516 


0.02507 


0.02451 


0.02447 



TABLE IV: Comparison of exact vàlues of the first 7 eigen- 
values of the CTOA operator for 7 = 0, 7 = ^ and 7 = f 
with the eigenvalues of the CTOA operator for 7 = 0, 7 = ^ 
and 7 = f computed using the Nystrom method. 
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